function Fig2_e()
lb = 25.66 / 6^0.5;
EC = 56.1 * 6^0.5;
scale_factor = 1/2/pi/lb^2;

% This proccessed data is processed by 
%  is located at .../processed data
result = readmatrix('processed data/result_one_half.csv');
r = result(:,2) * lb;
phi_r_total = result(:,3) * EC;
phi_r_imp = result(:,4) * EC;
phi_r_2DEG = result(:,5) * EC;
n_r_2DEG = result(:,6) / lb^2;


Fcolor='k';
Bcolor='w';
Fsize=18;  %font size

g = figure();
g.Color = Bcolor;
box on

yyaxis left
plot(r,n_r_2DEG)

ax = gca;
ax.Color=Bcolor;
ax.XColor=Fcolor;
ax.YColor=Fcolor;
xl=xlabel('r (nm)');
yl=ylabel('charge (nm-2)');
ax.FontSize=Fsize;
yl.FontSize=Fsize;
xl.FontSize=Fsize;
limit=[0,125,-1e-4,7e-4];
axis(limit)

yyaxis right
plot(r,n_r_2DEG/scale_factor)

ax = gca;
ax.Color=Bcolor;
ax.YColor=Fcolor;
yl=ylabel('nu');
ax.FontSize=Fsize;
yl.FontSize=Fsize;
xl.FontSize=Fsize;
axis(limit ./ [1 1 scale_factor scale_factor])
end
